This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 

BEST AVAILABLE IMAGES 



Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not hmited to): 



BLACK BORDERS 

TEXT CUT OFF AT TOP, BOTTOM OR SIDES 

FADED TEXT 
ILLEGIBLE TEXT 
SKEWED/SLANTED IMAGES 

COLORED PHOTOS 

BLACK OR VERY BLACK AND WHITE DARK PHOTOS 
GRAY SCALE DOCUMENTS 



IMAGES ARE BEST AVAILABLE COPY. 

As rescanning documents will not correct images, 
please do not report the images to tiie 
Image Problem Mailbox. 





WORLD INTCLLECTUAL PROPERTY ORGANfZATTON 

Iniernaiionat Bureau 




PCX 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) Intcrnutional Patent Clo^^Jlcation ^ : 

A61B 5/08, 5/05 



Al 



(11) Internalional PublicAlion Number: WO 97/37591 

(43) InlernntinnQ] PubHcoUon Date: 16 October 1997 (16.10.97) 



(21) IntcrnHlional Applicution Number: FCT/SG97/00013 

(22) Intcrnutional Filing Dalfc: 7 April 1997 (07.04.97) 



(30) Priurity Data: 
629,420 



8 April 1996 (08.04.96) 



US 



(71) Applicant (for all designated States except US}: RHEO- 

GRAPMTC PTE LTD fSG/SGJ: 17 Jurong Pore Road. 
Singapore 6] 9092 (SO). 

(72) Inventors; and 

(75) Inventors/AppIicDnls (for US only): SCHOOKIN. Sergei. 1. 
PIU/RUJ; Yurrovski por., 22 bid. 1. app 21. Moscow, 
11)020 (RU), 2UBENFC0. Viatchcslav C. fRU/RUJ; Iz- 
mailov,<:coic sh., 45, app 14. Moscow, 105187 (RU). BE- 
LIaEV, Konsttmtin R. (RU/RUJ; Foscivalnaya st. 5. app 
115. Moscow, 125565 (RU). M0R020V, Alckiandr A. 
[RU/RUJ; DezhnevasL 6 bid. 2, app 127. Moscow. 129642 
(RU). YONG, Wen. H. [SG/SGl; 43 Poh Htiat Drive. Sin- 
gapore 546628 (SG). 

(74) Aftent: CHAN & RAVINDRAN (PATENTS) PHE. LTD.; P.O. 
Box 2988. Raffles City. Singapore 91 17 (SG). 



(81) Designated States: AL\ CA, CN. JP. KR, RU. US. European 
paieni (AT, BE, CH. DE. DK. ES. FI. FR. GB, GR, IE, iT. 
LU. MC. NL. PT, SE), 



PubJishud 

WitJi iniernational search report. 

Be/ore the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
aniendrnenis. 



(54) Title: NON-INVASIVE MONITORING OF HEMODYNAMIC PARAMETERS USING IMPEDANCE CARDIOGRAPHY 
(57) Abstract 

A rneihod and appftrarus for dctcnninatioj\ 
of hcan race, heaa stroke volume, and cardiac out- 
put from thoracic bioimpedhnce signals and clec- 
trocordio^rams. A unique bioimpcdance electrode 
flrrungcmcnt is employed, and the bioimpcdance 
signals arc corrected for gain phase-frequency dis- 
loriion through the use of sinusoidal test signals 
through the measuring or detection electrodes [o 
identify distonions and correct for saroe during ac- 
tual measurements. Time-derivative bioimpedance 
signals are employed, the power spectrum cal- 
culated, and a novel autoconvoluclon procedure 
used to emphasize the heart rate hamionic. Breath 
waves and other signals not indicative of ilie pa- 
tient's cardiocycles are removed. Left ventricular 
ejection time is derived from the bioimpcdance sig- 
nals, and an improved version of Kubicek*s equa- 
tion is employed to derive heart stroke volume and 
thus cardiac output. 




FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to idencify States party to the PCX on the from pages of pamphlets publishing Inicmational applications under the PCT. 





Codes used to idencify States 


Al. 


Albania 


Its 


AM 


Anneniit 


Ki 


AT 


AuairU 


FR 


AU 


Au3irali;i 


ga 


A2 


Azerbflljnn 


CD 


UA 


Boiiila jrid Herzegovina 


C£ 


no 


&arbndoE 


GH 


HE 


fielgium 


CN 


DF 


Burkina Faso 


CM 


BC 


Bulgaria 


IIU 


HJ 


Benin 


IE 


DR 


Qrazil 


!L 


BY 


BcLsra^ 


IS 


Ca 


Cartada 


IT 


CF 


Central ArHcBn Republic 


JP 


CO 


Congu 


K£ 


Ol 




KG 


CI 


Cfric <J' Noire 


KT 


CM 


Cameroon 




CN 


China 


KR 


CU 


Cuba 


KZ 


C2 


CT/ceh Republic 


LC 


DG 


Ccrriitny 


LI 


UK 




LK 


Mr. 




LR 



Spain 

Finlnnd 

France 

Ciibon 

United Kingdom 

Ccorgio 

Gh^a 

GuioeB 

Greece 

nungory 

Ireland 

Israel 

fcclnnd 

linly 

Ken/A 
Kyrfyzeifiit 
DemocrAiic Peoplc'i 
Republic or Knren 
Rrpublic of Korea 

Sjinc Lucia 
Licchlenstcin 
Sri L;(nkii 
Liberia 



LS 

I.T 

l-U 

LV 

MC 

MO 

MC 

MK 

ML 

MN 
MR 

MW 

MX 

NE 

NL 

NO 

N2 

PL 

PT 

RO 

RU 

i»D 

SE 

SG 



Lcsoiha 


Sf 


Sloven in 


Lirhunnin 


SK 


2>lovnlcia 


Luxernbovirg 


SN 


Senejrol 


Laivio 


sz 


SwAzllutd 


Monneo 


TD 


Chad 


Republic of Moldova 


TG 


Togo 




TJ 


Tijikiiton 


The foirncr Yugoilftv 


TM 


'J'lirkmcfiiscoji 


Republic of M;u;cdonia 


TR 


Turkey 


Mali 


TT 


I'rinidad nnd Tahn^o 


Mongolin 


UA 


Ukraine 


Mouriinnia 


UC 


Uganda 


Malawi 


US 


Unlred Slates of America 


Men i CO 


L'2 


Uzbekiffrui 


Nipcr 


VN 


VIoi Nnm 


Ncihcflands 


YU 


Yugoslfivjft 


Norway 




Zimbabwe 



New Zenland 
Polind 
PonugAl 
Rcniiinlfl 

HuN^inn Fcdefiifioii 
Sudan 

Singapore 



wo 97/37591 PCT/SGJ)7/00013 

•1 - 

NON-INVASrVE MONITORING OF HEMODYNAMIC PARAMETERS USING 

IMPEDANCE CARDIOGRAPHY 



5 TECHNICAL HELD 

The present invention relates generaHy to cardiac monitoiing and specifically 
to the deteimination of hean rate (HR), heart stroke volume (SV), and cardiac 
output (CO) according to detection and complex analyses of thoracic bioimpedance 
and electrocaidiograph (ECG) signals, which permit precise detection of the start of 
10 left ventricular ejection. 

BACKGROUND ART 
Heart rate is the number of times the heart beats per minute. Hean stroke 
volume is the volume of blood pumped during each heart stroke. Cardiac output is 
the volume of blood pumped in one minute and is generally considered to be the 
15 most significant gauge of cardiac fimess. Physicians must frequently lely upon such 
cardiac parameters to diagnose heart disease, to zsscss a patient's overall health, to 
determine the most appropriate method of treatment, and to quickly discover sudden 
lapses in cardiac performance. 

The currently existing methods for measuring cardiac output and other 
20 cardiac parameters may be divided into two categories: invasive and noninvasive. 
The invasive methods require that a medical practitioner insert a measuring device 
into the patient's body, such as a catheter in the throat, and present numerous 
disadvantages to both patient and physician. The patient must often endure 
substantial pain and discomfort and the physician muse perform a relatively 
25 complicated procedure and occasionally expose himself or herself to the risk of 

contact with infectious blood. The noninvasive methods currently in use represent a 
major advancement, but still have signifvcant shoncomings. Most take 
measurements using ultrasound, phonocardiography, or electrical bioimpedance in 
order to calculate cardiac parameters. 
30 The methods which employ bioimpedance measurement involve placing a 

plurality of electrodes on a patient's skin (predominantly in the thoracic region), 
generating a high frequency, low amplitude eieciric current from cenain of the 
electrodes into the patient's body, measiirino tiic clianges in the electrical impedance 
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of the patient's tissue over time, and conrehiing the changw in cleciricaJ Impedance 
witJj cardiac parameters. 

The manner of arranging the electrodes on the patient's body plays an 
important part in the relative accuracy of the ultimate cardiac parameter 
5 tncasureiucms. Due to various ajiaroraical factors, electrodes must be placed over 
certain areas of the body to achieve optimum correlation between measured changes 
iii bioimpedancc aod cardiac paiaraeters. Many of the electrode configiinitions 
cunendy in use fail lo adequately take into account the paths foUowed by the lines 
of electrical potential thiuugh the thorax and thus create a distortion in the caidiac 
measurement. Moreover, a few electrode anrangements require the use of band 
electrodes, c.g.. influencing band electrodes A, B and measuring band electrodes C. 
D each having a width "n" (see HG. 1). These band electrodes typically wrap 
around a patient like a belt and further limit access to the patient, an especially 
undesirable condition during rcanimation procedures. The movements associated 
15 with respiration also make band electrodes very inconvenient when placed on the 
neck and chest. 

Perhaps the most significant problem with the presently existing 
bioimpedancc nietiiods is the imprecise raatheraaticaJ derivation of cardiac 
parameters from bioimpedancc measurements. The veniricular ejection time (VET) 
is a measurement of the time between the opening and closing of the aonic valves 
during the sysrole-diastole cycle of the heartbeat and it must be calculated as an 
intermediate step in determining cardiac stroke volume. The prior an does not 
teach a method for detennining ventricular ejection time with sufficient accuracy. 
Furtliennore, the prior an fails to account for the fact that VET is not a single 
25 event. In reality, there is actually a left VET and a right VET. Ii has been shown 
that the timc-dcnvative impedance signal is actually proportional to the peak aortic 
blood flow ejected by the left ventricle. The measuremenis of left VET and right 
VET for most pmients are generally very close, but even slight differences bet«/een 
them can create errors in bioimpedancc readings under the methods currently in use. 
30 Furthermore, the chssic algorithm for ejection starr time is elaborate, and 

works well only for healthy patients at rest. It is not accurate for patients under 
physical training or other physical .Lms.s, or for critically ill patients, such as those 
typically in intensive care units. 
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The conventional equation for deriving stroke volume from bioimpedance 
signals has become known as the Kubicek equation and is given as follows: 



SV = RiLjZ^^ • hZ 



25 



■ 

wheae SV is heart stroke volume, R is blood rcsistivity, L is the distance between 
the inner and outer voltage sensing electrodes, Zo is the mean tlioracic impcda/ice 
determined from the inner voltage sensing electnxles, and AZ is the impedance 
change due to blood influx. Kubicck's estimation of this value is 



10 where VET is the combined left and right ventricular ejection time, and (dZ/dt)^, is 
the maximum negative slope change of the time-diffecemiated impedance signal. 
Most bioimped;uice cardiac monitoring systems use some form of the Kubicek 
equation. 

Without further refinement, however, the Kubicek equation frequently gives 
15 inaccurate measurements. This is due in part to the fact that both ventricles 

contribute to impedance changes, and so Kubicek's calculated ejection lime (VET) 
cannot be associated with a particular, specifically the left, critical ventricle. 
Concurrently, Kubicek's AZ estimation becomes invalid when strong left-right 
ventricles asynchrontsm is observed. As a result, Kubicek's SV caJcuIaUon is often 
20 proportional to, but not equal to, the acmal heart stroke volume and must therefore 
be multipUed by some correlating coustanL. In addition, tlie prior art does not 
disclose a method for adjusting R in accordance with the fluctuation of a patient's 
hematocrit (red blood cell count). The adjustment of R i.s especiaUy important in 
patients undergoing blood infusion. 

Many of the mctlwds for bioimpedance cardiography require that the patient 
hold his or her breath during each measurement because respiration cjujscs 
interference in (he bioimpedance signal. Such methods arc inconvenient for some 
patiema and completely useless for other patients who are unconscious or otherwise 
unable to hold their breath. Some of the more recent methods include signal 
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processing c.ip..|biluy lo enhance the sigt,al. to identify the effects of respiration, and 
to eliminate defective signals so that errors are not introduced into the final 
calculations. Effective signal processing is gencmlly the key to insuring accuracy in 

bioin^P^nncc ctrirdiography and impfOvement5_in_this are^ cnn-r^pr^t^iem^eaii^ 

5 advances in the art. 
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DISCLOSURE OF iNVENnON 
The invention discloses a method of measuring hemodynamic paiametars 
using a novel combination of bioimpedance cardiography and electrocardiography 
iO Which allows medical practitioners to obtain an accurate, substantially continuous 
assessment of a patient's cardiac perfoimance. n,e bioimpedance and 
electrocardiogram signals arc measured over a common time interval of interest, 
preferably more than ten (10) heartbeats. 

The apparatus of the invention uses a series of spot electrodes adapted for 
placement on the surface of the patient's skin to generate a high frequency, low 
ampUtude electrical current through the thorax of the patient and to measure 
changes in the bioimpedance. The disclosed electrode configuration takes advantage 
of the physiological arrangement of electrical potential power fines in the body. 

Concurrently, a method of the invention may utilize elcctit)caitIiography to 
enhance, the accuracy of the ejccUon time detection from the measured bioimpedance 
signal. TJie cleclrocairiiogram can be obtained in any standard position, well known 
in electrocardiography. To eliminate the overall quantity of electrodes at the 
patient's sidn, this invention may use the same set of electrodes for bioimpedance 
and ECG measurements (see FIG. 6B). 

The invention may aUo involve the step of continuously adjusting the 
calculaaon of cardiac parameters according to changes in the rod blood cell count 
and according to vailing bodily composiitons of different patients, and can thus be 
used under a variety of different circumstances without a loss of accuracy. 

In addition, the ii.vention comprises a method of improved bioimpedance 
signal processing. It may employ a computer system to analyze the both the 
bioimpedance signal and the clectrocardiograjn in a variety of ways to pmvide an 
accurate repon of cardiac parameten. The computer .system may be used to make 
connections m the gain-pl,a.se-frequcncy characteri.stics caused by the transducers 
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u.sed to measure the bioimpetlance and ECG. The computer system communicates a 
sinusoidal test cuive to the transducer receiver and may then measure and record the 
sain-plmsc-frequcncy distortions created by the receiver. The "real" bioimpedancc 
signals subseqtiently received by the computer system throug h that fT:i n..riM r^.r mny 



paw tlirough filters which remove the transducer's characteristic disionions as weU 
35 brcaUi and movement aitifacts. Accuracy of gain-frequency response and phase- 
frequency response may he corrected to n/ithin 5 % . 

The invention may fiirther comprise a method of signal processing the ECG 
signal to detcimioe QR5 complexes (characterisuc heartbeat waveforms) and check 
10 point positions for use in refining die bioiinpedance signal calculations. The signal 
processing may involve sampling the measured ECG signal and filtering it to 
highlight the positions of the QRS complexes. Peak-to-pcalc amplitudes may be 
recorded and a threshold amplitude may be calculated. Dcfect-ftras QRS complexes 
are selected using the threshold amplihide. Finally, additional analysis of selected 
15 events may be jierfoimed to determine check point positions. 

The computer system of the inveation may derive and save in memory a 
time-derivative bioimpedance signal. The curve generated by the time-derivative 
bioimpcdance signal may be plotted with respect to time and represents repeated 
cardiocycles. The power spectatm of the bioimpedance signal may be calculated 
20 with discrete Fourier transforms and studied to estimate the patient's heart rate and 
to identify the fronts of each cardiocycle. A novel mathematical autoconvohuion 
procedure may be used lo emphasize the heart rate harmonic in the time-derivative 
bioimpedance signal. 

Breath waves may be removed by generating an envelope within the power 
spectnjm in which cardiocycle signals should be found and then removing those 
signals which He outside the envelope. 

The computer system may employ a new method of deriving the effective 
left ventricular ejection time (ELVET) from bioimpedance signals and ECG signals. 
In particular, precise detection of the ejection stan time is based on complex 
analyses of the bioimpedance and ECG signals. The derivation of ELVet is a 
detailed, multi-step analysis which ijivolvcs finding various points on the time- 
derivalivc bioimpedance curve based o.i corresponding points on the ECG curve and 
detennining whidi of these pomts mosi accurately reHect cardiac events. It rcqi,ir«.s 
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making a variety of maihcmatical calculations, including .nalung approximacioiis for 
differentials ol" the time-derivative bioimpedance curve. 

In heallhy patients there is very Uttie left-right ventricular asynchit)nism or 

^synchron ism betw een opening of left and ri ght ver^tririP. valves. However, in 

5 seriously ill patients, tlie Kubicek equation generally underestimates the time- 
derivative bioimpedance signal due to Jeft-right ventricular asynchronism. The 
computer system, in accordance with the invcncion, automatically adjusts for signals 
which exhibit asynchronism between opening of left and right ventricle valves to 
calculate the ccrrect value ofAZ, the impedance changes due to blood influx. 

After the computer system calculates ELVET ajid AZ. it employs an 
improved fom of Kubicck's equation to derive Heart Stroke Volume from ELVET. 
blood resistivity, die patient's specific body constitution, and the maximum 
bioimpedance change. The computer system analyzes the time-derivative 
bioimpedance signal to eliminate cardiocydes with certain aberrations. 

Finally, the invention fvrtlier comprises a method of detecting valid 
cordiocycles. 



25 



BRIEF DESCRIPnON OF DRAWINGS 
FIG. I i.s a schematic depiction of a band-type electrode arrangement on a 

20 patient; 

FIGS. 2A, 2B, and 2C schematically depict spot-type electrode placements 
according to the present invention; 

FIG. 3A depicts anatomical zones of interest for placement of measuring 
chest electrodes according to the present invention; 

FIG. 3B illustrates the posiUoning of the anatomical zones. of FIG. 3 on a 
patient; 

FIG. A illustrates the stiucturc of an individual chest electrode according to 
the present invention; 

FIG. 5A IS a graph of the gain versus frequency characteristics of an A-fllter 
30 according to the invention employed for signal differentiations and harmonics 
suppression; 

FIG. 5B is a graph of the gain versus frequency response of a B-filier as 
employed in tl,c presem invention to highlight the froiit.s of cardiocydes; 
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RG. 5C: is a graph of filtration employed to rcmovd breath waves' power 
spectrum from the summary (breath and cardiosvave) power spectniin in the presam 
invention ; 

P^"G- 6a is a brancli chart of the methodol og y of the present inventio n: 

5 FIG. 6B is a flow chan of signal processing and depicting the apparatus 

employed in performing the method of the present invention; 

FIG. 7A depicts tbe gain-frequency characteristic G(f) of a bioimpedance 
transducer according to the present invention; 

FIG. 7H depicts the phase frequency characteristic P(f) of a bioimpedancc 
10 transducer of the present invention; 

FIG. 7C depicts the gain frequency characteristic ^ of an "R" or 

restoring filter corresponding to the characteristic of FIG. 7A as used for GPF 
collections according to the instant invention; 

FIG. 7D depicts the phase frequency characteristic -P(f) of an "R" or 
15 restoring filter corresponding to the characteristic of FIG. 7B as used for GPF 
corrections according to the instant invention; 

FIG. 8A depicts the resulting gain frequency characteristic G(f) of a 
transducer signal as corrected by an R-filter according to tho present invention; 
FIG. 8B depicts the resulting phase frequency characteristic P(f) of a 
20 transducer signal as corrected by an R-fiJier according to die present invention; 

PIG. 9 is an illustration of a time-derivative bioimpedance signal generated 
according to the present invention, plotted with respect to lime and depicting 
exemplary impedance changes due to blood influx (AZ) over the same time period; 

FIG. 10 is a graph of the electric potential traversing the hean of a patient 
plotted with respect to time (electrocardiogram), depicting the peak-to-peak 
amplitudes (E, - front, E, - back) of the QRS complex in a single heanbeat; 

HG. 11 is a scatter diagram of all peaks measured in an electrocardiogram 
over a period of ten seconds The coordinates of the graph arc (E,, Ej), ami the 
QRS complexes are within ihe doued circle in the lower right Jiand comer; 

^G- 12 is a graph of the ECG signal and the time-diffcrcntiatcd 
bioimpedance signal plotted against time ijidicaling check-points P. S.,. and S, on the 
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ECG Signal, relative to check-points Q, S, and A on the time.-dirferentiatcd 
bioimpedoncc signal; 

HG. 13 is a graph of ECG, bioinipedance Zft), and time-diffeientiatcd 
bioimpedancc dZ/dt, signal.': with respect to time indicating the bioimpedance 
difference between the S and Q points. in a healthy patient; 

FIG. 14 is a graph of ECG. bioimpedance 2(0, and time-differentiated 
bioimpndancc dZ/dt, signals with respect to time indicating the bioimpedance 
difference between the S and Q points, Z^, in a patient with ischemic hcajrt disease; 
and 

FIG, 15 is a graphic depiction of the correlation of the inventive method of 
detennining ELVET to an idea], ultrasound determination in comparison to the 
classic Kiibicek algorithm methodology. 

BEST MODES FOR CABRYING OUT THE INVE^mON 
The first step in the present invention involves taking bioimpedance 
measurements over segments of tissue on a patient's body. Electrodes must be 
placed at appropriate points on the surfece of the skin to generate a high frequency, 
low amplitude electric current and to detect changes in the generated canent after It 
passes through the segments of tissue (see HGS. 2A, 2B, 3 and 4). The electrodes 
arc "spot electrodes" rather than "band electrodes" in order to maximize the free 
area on the patient's body. The spot electrodes are preferably of the disposable, 
one-use type. The patient thus has increased freedom of movement and medical 
practitioners have more access to the patient's skin for other medical procedures, 
such as tlie introduction of catheters and the administration of anesthesia. 

The bioimpedance electrode system employs a total of six electrodes: a pair 
of detecting (measuring) electrodes 20 at the xiphoid process level, a pair of 
detecting (measuring) electrodes 22 positioned lateraUy on the neck, an influencing 
electrode 24 on the left leg, and an influencing electrode 26 on the forehead. 

The influencing electrodes 24, 26 may be standard ECG spot electrodes with 
a contact area of 2 centimeter.'; by 2 centimeters (cm). The upper influencing 
electrode 26 is preferably placed on the middle of the forehead, at the mid-line 
thereof. The lower influencing electrode 24 is prefarably placed on the left knee or 
somewhere below the left knee such as (he lefl foot. If necessary, the lower 
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influcncing elecuacle 24 may also be placed above the knee level, provided thai the 
following condiLion is satisfied: 



L> SJt, 



where L Is Uie distance between influenciiigr electrodes and R is the radius of the 
5 chest. Tl)e left leg is used instead of the right leg to account for the anatomic 
asymmetry of die hean. The physiological positioning of the aortic aich. through 
which a significant concentration of electrical potential power lines pass, makes the 
left leg most suitable for the lower influencing electrode. The arrangement of 
influencing electrodes in this manner guarantees the uniform distribution of 
10 influencing cunent power lines between measuring electrodes and thus helps to 
mim'mize the cxror in the final cardiac parameter measurements. 

The upper pair of measuring electrodes may also be standard ECG spot 
electrodes with a contact area of 2 cm by 2 cm. These electrodes are placed 
symmetrically Jilong the lateral lines of the patient's neck about the perimeier of the 
15 patient's neck 27, a distance S above the base of the neck. The distance S is 

defined as the distance between the base of the neck and the cenier 302 of electrode 
22 and is preferably approximately 4 cm. The base of the neck is defined to be 
located at the point of maximuni curvature of the lateral lines of the neck. Placing 
the upper measuring electrodes 22 in this area avoids the error that would otherwise 
20 result from the nonlinearity of electrical power lines at the neck-chest junction. 

The lower chest pair of measuring electrodes 20 each have a contact area of 
12 cm^ to 30 cm^ If this contact area is either reduced or enlarged, the heart 
stroke volume will be underestimated A contact area of less than 12 cm^ provides 
insufficient depth of measurement, a panicularly serious problem with larger 
25 patients, and a contact area of more than 30 cm^ causes the measuttiment to extend 
into additional anatomical regions. 

The individual chest electrodes 20 are each preferably comprised of a set of 
four standard ECG spot electrodes 28, each with a contact area 29 of 2 cm by 2 cm 
with the top pair 34 of spot electrodes 28 of each chest electrode ai xiphoid process 
30 level 38 (see FIGS. 2B and 4). All contact areas 29 are connected with foil or wire 
31. Tlie distance G separating adjacent spot electrodes is approximately 5 cm. The 
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concact areas 29 are placed on the body using a conductive gcl. if noi an iiiiegrai 
feamrc of the .spot electrodes 28. Such a design takes advantage of Uie anatomic 
features of zones I and HI of tfie body (see HQS. 3A and 3B) and insuras a 
minitnum of eiror since the measurement is taken at m adequate depth and allows 
5 for variations in the bodUy constitutions of different patients. The chest dcctiodcs 
20 as describctl are placed laterally on opposite sides of the chest (sec FIGS. 2A, 
2B, and 2C) at the xiphoid process level 38. As shown in FIGS. 3A and 3B. 
placement is ijt zones H and IV, within an area anteriorly or posteriorly plus or 
minus about 10 centimeters of a lateral line extending tluough the body at xiphoid 
10 process level 38. 

The infltienciug electrodes 24, 26 generate a high frequency, Ioav amplitude 
cuttent into tho paticnfs body and the detecting electrodes 20, 22 measure the 
current after it passes through body tissue. The electrical impedance of the tissue 
can readily be determined from the difference between tlie generated current and the 
measured current. The electrical impedance of the tissue varies over time as a 
result of blood flow, respiration, and other factors. 

The present invention aho uses ECG signals concurrendy measured with tlic 
bioimpedance signals. With the exception of post-measurement signal processing to 
remove hardware artifects described below, the ECG signal measurement is 
perfomied in a conventional manner, and so will not be further described. 

A preliminary step in the present invention is to determine the f ain-phase- 
ftequcncy (GPF) characteristics of the analog input devices (ECG converter and 
impedance converter) for subsequency use in signal processing. 

A computer systejn. which may comprise a specifically programmed general 
purpose computer such as a personal computer, receives the electrical current 
measurements from the detecting electrodes 20, 22, determines both the impedance 
of the intcrcedi/ig tissue (bioimpedance) and the electric potential traversing the 
heart (ECG) as n function of time, and uliimaiely calculates HR (heart rate); SV 
(heart stroke volume), and CO (Card.ac Output). The equation for determining SV 
30 is given as follows; 
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SV= K 'f '(LI Z^^AZ 



where K is a nnvel scale factor, P is Uie specific bleed resistivity, L is the distan.cc 
between the voltage measuring or sensing electrodes 20 and 22. is the mean or 
base thoracic impedance (determined from sensing electrodes 20 and 22) and AZ, 
the impedance changes due to blood influx. AZ is calculated as 



AZ = ELVEr*(dZ/di)^^Z^^, 



where ELVET is the effective left ventricular ejection time, (dZfdt)^ is the 
maximum absolute value of the tijne-differentiated impedance signal obtained from 
the two measuring electrodes, and Z,, is a novel coirsction factor that lakes into 
account left-right ventricles asynchronism, equal to the bioimpcdance difference 
between S and Q points (see HGS. 13 and 14 and accompanying text). 

The K tcmj accounts for variations in body compositions of different 
patients. Li order to obtain a value for the IC feoor. the medical practitioner first 
measures the height and weight of the patient and the perimeter of the patient's neck 
27 and chest 36 using conventional means or coramercially available ultrasonic 
measuring means. The medical practitioner then inputs these values into the 
computer system, which in tum uses the values to compute the effective cross- 
sectional area of the chest and the K factor. The effective cross-sectional area of 
the chest (SCHES17 is given by: 

SCHEST = {PCHESP- * {PNECK • PCHEST) * PNECK^-) / 12n, 

where PCHEST is the perimeter of the patient's chest 36 and PNECK i% the 
perimeter of the patient's neck 27. Then the K factor is calculated as follows: 
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wlierc 11 is thi: patient's height, W is the paUent's weight, and ICo, Kj, K^, K3 are 
gender and age dependent and lie in ranges" 

J^o ^ [1-4]; A'^ <• [3-16]; e [0-1]; e [0.1-2] 



The present invention can thu.s be used on patients of varying body consiicutions 
5 without a loss of accutacy. 

Tlie elec:trical resistivity of human blood is not a constant. It varies among 
different individuals and even in the same individual at different times. Blood 
resistivity is paiiicularly susceptible 10 flucluaiion in patients undergoing blood 
infusion. As a result, an accurate system of bioimpedance cardiography must 
10 include means Jbr continuously modifying the blood resistivity tejm of the Kubicek 
equation. 

A patient's specific blood resistivity depends largely upon his or her 
hematocrit. The relatiouship between these two values for capillary blood is as 
follows : 



15 



20 



where P is the specific blood resistivity and H is the capillary hemaiocrit. This 
relationship is adopted from V,L Arinchin el aL, "Taking into account electrical 
resistance of blood will increase accuracy of chest tetrapolar rheogiaphy method, " 
Journal of Pediatrics (U.S.S.R.) 1987, v. 7, ppJ9-52. Tlie hematocrit can be 
measured using any commiircially available method. It may either be inputted into 
the computer system by the medical technician, or sent by electronic means directly 
from the hematocrit measuring device. 

The invention utilizes a novel method for processing the ECG signal, (after 
hardware artifact removal) comprised of the foUowiiig steps: 
0) signal approximation from sampling points, 

(ii) special filtering co highlight the positions of QRS complexes, 

(iii) measuring peak-io-peak amplitudes for a given, recorded dme 
interval, 

(iv) calculation of amplitude threshold, 
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(V) QRS selcctioa wixh the calculated amplitude ilireshold, 
(vi) acldiiioiial analysis of selected events to deteraiinc chedc point 
positions. 

Step (i) al>ove is desirable to increase the accuzacy and reliability of QRS 
5 detennination unless sampling frequency is extremely high. As a matter of 

practicability, such high frequency sampling is undesirable as consuming excessive 
processing time and memory capacity, and is impractical to effectuate the signals of 
restricted power spectrum. It is contemplated that two apprDache.s lo ajpproxiinacion 
may be suitable for use with the invention. First, it is well known that any signal 
10 s(t) with finite spectrum (and defining the highest haimonic as oi^^lirfj, is ftiUy 
described by its samples at points s(n*AT). where AT<l/2f„ is the sampling 
period, and n is an integer. A precise approximation of such a signal is given by 
the equation: 

sic) = Y, sinAl) /^^, ^ = E ■ 



'(at"''*, 



where <p^(0^j/ni:j^Tr- 

The same result can be obtained in another way by first caJculating a Fourier 
transform of the signal 5(t) and adding a small phase shift Ay>, lo all hairaonics, so 
that 

A(p/2t^ = t = canst, V^. 

Then, after calculation of the inverse Fourier transfonn, ihe approximated values of 
20 s(t) at the points shifted by ^ from the original samples arc obtained. This latter 
approach is more efficient in calculation. Both methods permit lower sajnpling 
rates that lead to less consumption of memory and can be employed on demand to 
calculate a precise approximation of the original signal. 
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Tlic next step is filtering ihii ECG signal lo highJigJil the positions of QRS 
complexes. FIG. 10 illustrates a typical QRS complex, a signal peak with greatest 
amplitude measured from peaJc-to-peak in a single heartbeat. A symmetrical fmiie 
impulse response (PER) digifil filter is calculated from the desired gain-frequency 
5 characteristic (TGCfi,^). The desired GFCnu^ is elaborated from analysis of power 
spectrum of QRS complexes and has a passband from 6 Hz to 22 Hz with the 
maximum at 12.5 Hz. Using a discrete Fourier transform, the desired GFC^hv is 
conveitixl to fmite impulse characteristic according to an algorithm for filter 
synthesis described in V.S. Gumikov, "Filtration of measured signals/' Leningrad. 
10 Energoatonuzdat (USSR) 1990, pp. 172-181, incorporated herein by reference. 

This filter passes the QRS complexes and suppresses breatli and movement artifacts 
in the ECG signal, also as P and T waves. 

The next step is to calculate the peak-co-peak amplitude threshold and select 
valid QRS complexes. The computer system measures each local peak of the 
15 filtered ECG signal by its front (E,) and back (E,) amplitude fronts, sec FIG. 10. 
For each local pealc, E, is the measured from the peak's anterior, or leading local 
minimum to the next nearest maximum, and Ej is measured from the peak's 
maximum to its posterior, or trailing, minimum. FIG. 11 depicts a distribution 
(scatter diagram) of peaics by their (E^, Ej) coordinates for a time interval or period 
20 of 10 seiconds. FIG. 11 also show5 the QRS complexes highlighted within a dotted 
circle in the lower right hand comer of the figure. Each pealc is chaiaclcrizcd by 
its (Ej, E2) vector arid amplitude Aj, where 



The computer system Ihen searches the sorted {A^} array for the maximum 
difference between A, and A^^,. If liic maximum is found for an exemplary k-th 

25 element, then (he amplitude threshold is calculated as J = (Ai^+Aj,^j)/2. Thus, the 
QRS complex is detected in point j if A^ exceeds the threshold T. The threshold T 
is adapted for each 10 second block or incen/al of ECG data as T. = T.^.j 4- aT^, 
where T^^., is the adapted threshold for the previous data block. T^ is the calculated 
threshold for the current block, and a is a parameter of adiipution in U\e range of 

30 0-1. 
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Eich QRS cojHpIcx identified by the use of ihc above "threshold" 
methodology is ihen further aiwlyted within an interval of -50 to +200 milUseconds 
(ms) from the dciermined QRS position. For each QRS complex identified, the 
computer system determine* Uie amplitadc. sequence of peato. and derivative of 
peak fronts to arrange 3 clieclc points (HG. 12): the start of the QRS complex (P- 
point), tlie maximum deviation from thp. base-line CS»-point. which coincides with 
peak R in a nonnal ECG), aud the end of the QRS complex CS,-point). Tliese 
check points arc used to refijic the analysis of the bioimpcdance signal described 
below. 

The invention utilizes a novel method for processing the faioimpedance signal 
comprised of the following steps; 

(i) tligital filtration and phase correction, 

(ii) heart rate estimation, 

(iii) suppression of breath waves, 
^5 (iv) determination of cardiocycles, 

(v) arrangement of check points, and 

(vi) selection of cycles without interference artifacts. 

Hie first part of the electronic filtration involves passing the signal through a 
"restoring" R-filter to achieve gain-phase-frequency (GPF) correction. The R-filter 
compensates for the distortions caused by the particular electmnic transducer that is 
used to measure bioirapedance changes. It is wcD known that GPF characteristics 
of a bioimpedaace Uansducer (KG. 7) may grcaUy influence the shape of the 
bioimpedance curve. These interferences must be removed from the signal. The R- 
fllter uses posterior signal processing to connect hnear GPF distortions. It is 
constructed in such a way that the system of the bioimpedance transducer plus R- 
filter has GPF characteristics with zero phase shift and constant gain at the given 
range of frequencies, for example from 0.3 Hz to 30 Hz for the bioimpedance 
signal (see FIG. 8). Tlius. correlation of tlie outputs of different bioimpedance 
devices may be achieved. With the exception of the filter parameters and opeiati;,g 
chancteristics specified as desirable or critical for the R-filtcr and other filters 
described herein, constnjction of same is conventional and within the ability of 
those skilled in the art. and so will not be further described. 
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nic first step of the CPF correction mvolvcs connecting the bioimpedance 
transducer with a source of an elecirically-generat«yj sinusoidal impedance signal 
and then measiiring the output' from the transducer. The clcctronicaHy-gencrated 
sinusoidal impwlance signal has an amplitude of 0. 1 Ohm to 0.2 Ohm ^vith respect 
5 to a baseline, for example, 100 Ohm to 20O Ohm. Such a sigiial has been 

developed using a voltage-to-impedance converter consisting of a photoresiaior, a 
photoeniitter (photo diode), a pawer supply, and an analog-digital-analog (ABA) 
computer interface. The paired photoresistor and photoemitter are coupled inside a 
light-protected housing so that the photoresistor changes its impedance according to 
the photoemitter's light intensity. The ADA conversion process includes digital-to- 
analog conversion of the mathematically-modeled sinusoid with a frequency of 19 
IcHz and aualog-to-digital conversion with a frequency of 100 Hz. with 12 bit 
resolution. Through the interface, the computer produces a set of test sinusoidal 
signals with frequencies in the range from 0 Hz to 75 Hz and records the responses 
of the transducer. The operating charactoiistics of the voltage-to-impedance 
converter include an input signal of 0 V to 5 V, an output signal of -0.1 Ohm to O.I 
Ohm with a baseUne of 100 Ohm to 200 Ohm (as previously noted). The qPF 
characteristic H(f) of the transducer may then be calculated from the spectnims of 
the initial test signals and the resulting responses of the transducer and presented as 
a graph or stored in an ASCH or other memory file. The system uses the calculated 
GPF characteristic H(f) of the transducer to calculate the "restoring" R-filter. The 
GPP characteristic of such R-fiiter can be formaUy written as 1/H(f) in a certain 
frequency range. The R-fiUer also provides frequency bounds through low and high 
frequency filters to provide suppression of random low and high frequency 
interference (see FIG. 7). The filtration ivith the R-filter may be done in the 
frequency domain using a Fourier tmnsform. It is preferable to use a Gaussian 
window with Fourier transfomi to eliminate the boundary effects of the recorded 
signal: 

Git) = exp[-2(or/(27))^] 

where 2T is a duration of recorded signal, t < T i.^ a time, and -a- is a predefined 
constant prcferal)ly in the range of 2.5-3. Multiplication of the Fourier image of the 
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rccordcrj sJgnal by the R-filtcr's GPP chiimcteristics rcsulcs in suppression of GPF 
distonions and additional filtration of the aignaJ. Reverse Fourier transformation 
and division l,y (he Gaussian window may also be employed. The same steps may 
be employed in the time domain without Fourier transforniation. The signal after 
5 R.fiitiation is referred to as the "restoied" signal, and this signal is used for further 
calculations. 

It should also be noted at this time that the identified GPF chaiacteristics of 
tlie ECG connection are processed to remove hardware artifacti from the ECG 
signal ill a similar manner to diat described above for the bioimpedance signal. 
GPF correction of both the ECG and bioimpedance signals promotes toie 
correspondence of time intervals and event times between the two signals. 

Tlie next step in the bioimpedance signal processing is heart rate (HR) 
estimation. The present invention uses two ways to calculate HR. The oidical way 
is to detect R-pealcs on the ECG signal as described above and calculate R-H 
interval. The inverse value multiplied by 60 corresponds to heart rate. If ECG 
signal cannot be processed to detect R-peaks for some reason, the second way is 
used. In the second way. the power spectmm of the "restored" bioimpedance signal 
is calculated with di^cnate Fourier transform and used to estimate the patient's heart 
rate (HR). Veiy oi^en a breath hannonic is the biggest one in the power spectniro 
of the bioimpedance signal. Consequently, it must be suppressed and the HR 
frequency response highlighted. A special tiansformation is nsed for this purpose. 
First, the power .spectrum (PS) of the "rtistored" signal is multiplied with the gain- 
frequency characteristic of the A-filter (see HG. 5A). This filter differentiates the 
signal and addiiionally suppresses harmonics below a certain frequency preferably 
selected at a rr-n^je from I Hz to 3 Hz because breath wave harmonics commonly he 
below 2 Hz and HR haniionics above 0. 8 Hz. The power spectrum of clear carriio 
signal consists of repeated peaks at frequencies HR, 2*HR, 3*HR, etc. 
Consequently, the following autoconvoluiion of the power .spectmm will emphasise 
the heart rate harmonic: 
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whcre ASl(i) is a result of ihc auioconvoludon of the power spectnmi and PSa(i) is 
Ihe power of a )!ivcn spectral line with frequency i which previously passed tJiough 
the A-fiUcr. It is prefeiable to have in the above product only the first tlirec 
elements PSa, btxausc the liighcr the frequency is, the lower the signal-to-noise 
5 ratio. The computer system uses the autoconvolution to search for Mas, the 
maximum value of ASl(i) in the nmge of 0.6 Hz to 5 Hz. The frequency 
associated with Mas is regarded as an estimation of HR. The estimation of HR is 
then used in ad<litional filtration (see HG. 5G) and cardiocyclc-identification 
procedures. 

10 The aberrations in the bioimpedance signal caused by respixauon must be 

removed to increase accuracy and to insure proper identification of the cardiocyclcs. 
Usually, breath frequency is less than heart rate frequency, but breath waves create 
a power spectmm which overlaps the lowest harmonics of the power spectrum 
created by the cardiocycles. So, it is impossible to remove the breath waves' power 
15 spectrum from the summary power spectmm entirely (see FIG. 5C). Cardiac 
strokes are a more slable,^ repetitive process in comparison with breath. 
Consequently, we can consider that their power spectrum consists of several narrow 
pealcs. All power spectmm harmonics between main cardiowaves' spectrum lines 
are combination of lateral slopes of these main spectrum lines and a noise power 
20 spectrum. If these internal harmonics are decreased, noise is mainly suppressed, 
and cardiowaves only slightly. As subsequently noted, these ideas form the basis 
for the algorithm for breath wave filtration. After estimation of the heart rare (HR), 
as described above, the first and second harmonics of the cardiowaves* specirum art: 
determined. The local minimums at the power spectrum nearest to these spectrum 
25 peaks can be considered as their bounds. All hannonics below the lower bound of 
the second peaJc except for those within the fu^st peak bounds are multiplied by a 
predetermined value less than 1 (as schematically shown in FIG. 5C). This results- 
in elimination of bnsaih wave amplitude, because the lattcr's power spectrum lies in 
the multiplied 7x^ne, but only slightly affects the cardiowaves. 
30 The next step in bioimpedance signal processing is cardiocycle detection. 

The invention also uses two ways (o do this. The fim way h to place the 
cardiocycle borders according to QRS complexes posidons in the ECG signal. If 
the ECG cannot be processed, the second way is used. In the second way, the 
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" restored" sifjnnl as menLioriud above is passed through the B-fiJc&r (FIG. 5B) to 
highlight the (nmis of the cardiocycles. This is a differentiating fiJter with sinusoid- 
like frequency bounds. The pass band of the B-filter is adapted for frequency 
harmonics that produce a main contribution to cardiocycle front (preferably frDm 1 
5 Hz to 6 Hz). Tlie area between a local minimum and the next local maximum at 
the signal passed through the B-fiher is regarded as a cycle front and described with 
a peaJc-lo-peak change in time and a pealc-to-peak change in amplitude. The 
computer system then generates a time-amplitude envelope by analyzing the fir^c 
five (5) lumnonics of the power spectrum generated by the signal after it passes 
10 through the B-fiJtcr. The cycle fronts are examined and certain ones are marked for 
further analysis if their peak-to-peak changes in time and peak-to-peak changes in 
amplitude are within the time-amplitude envelope. To increajse the reliability of 
caidiocycles recognition, the computer system calculates mean and variance of 
peak-to-pealc auiplitudea for die selected regions, 
1^ If the variance-per-mean ratio is less than a predetermined value, preferably 

0.3, then all of the marked fronts are transmitted to the next stage which involves 
the arranBemcni of check points and the sclcsction of defect free cardiocycles. In the 
alternative, if the variance-per-mean ratio is greater than a predetermined value, 
then additional analysis mast be performed. The additional analysis comprises the 
20 following steps: (i) the regions under examination arc separated into two groups 
according to their pcak-io-peak amplitude: those above the mean value ("the upper 
group") and those below the mean value ("the lower group"); (ii) the means, Ml 
and M2, and variances or standanl deviations, VI and V2, are calculated for each 
group; (iii) the appropriate values are inserted into the following inequality for each 
25 group: 



where a is a predetermined value, preferably 1.96, and V| and V2 are standard 
deviations for the lower and upper groups respectively; and (iv) if the inequality 
holds true, then the regions in the upper group are taken to be the fronts of the 
cardiocycles ami the regions in the lower group are eliminated from further 
consideration, otherwise all selected regions proceed to the next stage. 
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riie computer system identifies ccnain check points in tlie time-diffcrcntiatcd 
bioimpcdance signal to calculate the effective left ventricular ejection time,. ELVET, 
as a prcliminai7 step in determining hean stroke volume (SV). 

The present invention uses ELVET, a term which rcpicaente only the left 
5 ventricuhir ejection time, rather than Kubicek's VET, a tctm which nspresents the 
combined left and right ventricular ejection times. It is known that the value of the 
time-dilTorentintcd impedance signal is proportional to peak aortic blood flow 
ejected by the left ventricle. ThersfoTB, the most accurate calculation of SV 
requires that LVET be used. LVET is calculated from the following equation: 



10 



where ELVET is the time between the moment the left ventricular valve opens 
(S-point) and the moment thai it begins to close (T-point), and LVPT is the 
protodiastoly time (the time it takes for the left ventricular valve to close). - LVPT is 
not readily detectable with hemodynamic monitoring means because the changes tn 
15 blood flow are insignificant during the protodiastoly time. For this reason, the 

present invention uses the product of ELVET and (dZ/dl)ovix. the maximum absolute 
value of the first chest impedance derivative with respect to time, to account for the 
lack of an LVPT measurement. 

The calculation of ELVET requires an analysis of the curve generated by the 
20 graph of Y(x), time-derivative bioimpcdance, plotted with respect to lime (se^ PIG. 
9), The computer system first finds a global maximum of lime-derivative 
impedance, Y(x), over a given cardiocyclc and designates it as point A. The 
computer system then traces back in time from the A-point to the point in time 
corresponding to the S.-point on the ECG signal, and looks for abnormalities in the 
25 bioimpedance signal between those two points (see FIG. 12). The abnormalities of 
interest are: (1) dZ/dt zero crossing, (2) local minimum in dZ/dt, and (3) local 
maximum in the third derivative of the bioimpedance signal, cPZ/dt^, If there is no 
abnormality in liie bioimpedance signal found in the time interval between points A 
and S„ the cycle is considered defective and rejected from further consideration. If 
30 any abnormalities are foimd to the right of S,, the abnormality closest ro S„ 

approached from the right is selected as the ejection siaa lime, S. Otherwise, the 
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abnonnality closest to Sh approached from the left is selected as the ejection surt 
time, S. The usage of the ECO signal increases stability of the S recognition iji 
complex cases. 

To identify the end of ELVET, point T, the corapuier system fmt finds 
5 point To, which is taken to be either the Hrst or the second local minitnum after 
point A at tiie time differentiated bioimpcdance signal, labeled as T| and T} 
respectively (sec FIG. 9). Tlic computer system chooses between T, and T, after an 
anaJysis of die depth (amplitude) of the curve at each point. If dq3th of the second 
minimum is greater than a predetermined fractional value of the depth of the first 
10 minimum, then T^ is selected as Tq. Otherwise, T, is used. To increase the 

stability of Tq detection in noisy signals, the invention looks for the back ot trailing 
edge of die T-wave in the EGG signal. The back edge of the T-wave is detected 
from the local maximum ac EGG signal next to QRS complex up to the next local 
maximum at tlic graph of curvature of EGG signal versus time. If one of the points 

15 Tj or Tj is out of bounds of the T-wave' s back, the other point is used as Tq 

regardless of its amplitude. The "effective end" of the ventriculai' ejection, point T, 
is then identified as the nearest local minimum before point Tq on the graph of the 
curve generated by the second derivadvc of Y(x). ELVET is calculated as the tinie 
distance between points S and T. 

20 The invtintive method of determining ELVET has been correlated to an 

ultrasound ELVET determination with a correlation coefficient of r = 0,86. 
Kubicek's classic algorithm g\yn only a correlation of r = 0.71. See FIG. 15 for a 
graphic depiction of the correlations of die classic versus the new methodology to 
the ideal. Con.<;equently, and in contrast to Kubicek, left ventricular ejection dme as 

25 measured by the inventive method is measured in substantial isolation from the right 
ventricular ejection time. 

Normally, the point of the ejection start, S, coincides with the zero crossing 
of the time-differentiated (dZ/dt) bioimpedance signal (sec FIG. 13), but this is not 
the case for seriously Dl patients (see FIG. 14). For such patients the point of the 

30 ejection start is often placed at the "stair" or "abnormality" of the first front of die 
bioimpedance signal. This "stair" or "abnormality" characicristic (also lermed a 
"prewavc") of an ill paticnl is commonly referred to as left-righr ventricular 
asynchronism. 
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For heaJiliy patienLs, the impedance value at ihc beginning of QRS complex 
at ECG (Q-poiiU), is almo.st. the same as the impedance value at the S-point, Z,. 
This impedance, difference ix measured as Z^^, (^cc FIG. 13). However, for ill 
patients, Zp.^ can be significant, (see FIG. 14). Where Z^^ is small, Kuhdceic's 
5 estimation of dZ=(dZ/dt)„pjt*VET, is fairly accuraie. For seriously U) patients 
exhibiting left-right ventricubr asyndircnism, however, the Kubicck equation 
generally underestimates AZ. Thus, the premise of the "prewave" should l?e 
accounted for in maJcing the calculation. This invention compensates for the 
inherent underestimation of AZ in ill pacients by adding the bioimpedance difference 
10 between S and Q points, Z,.^ to Kubicck 's AZ estimation. By making this 
compensation, AZ is more accurately estimated, and the regression between 
bioimpedance and therraodilution cardiac output values is linearized, TIius, the 
computer system estimates AZ as 

A2 = (dZ / - ELVET ^ Z,,^ 

15 where Z^^ is the bioimpedance difference between S and Q points. 

After the computer system arranges aU check points, it eliminates 
caTdiocyclcs with certain aberrations. Fuzzy logic and fitness algorithms may be 
employed in this procedure. Several criteria are used for this purpose. The 
computer syscem first confirms that the time distances between the points describe<I 

20 above for each cardiocycle (e.g., point A, point T, Point S) do not exceed certain 
bounds. It also verifies that the amplitude difference between the start and die end 
of the cardiocycle do not exceed a predefined value- The amplitude change between 
the start and end of the cardiocycle must not exceed a predefined percentage of the 
maximum amplitude over that cardiocycle Furthemnore, the ratio of the time- 

25 derivative bioimpedance signal amplitude at point A i:* the time-derivative 

bioimpedance signal at poinr T muse be greater than a preiletermined value; that is: 
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where Y(A) and YfT) are the values of the time-derivative bioimpedance signal at 
points A and T respectively, and "c" is the predetermined value. All cardiocycIe.s 
that pass this stage are considered as "not very bad/' Tl)e computer system iben 
checks for the "neighbors criterion" to eliminate the effect of random noise in the 
5 bioimpedance signal. A three-dimensional "nearness" vector with elements (Ali, 
A2i, A3i) is calculated for each pair of cardiocydes in a 10 second time block. The 
individual elements of the nearness vector are dctennined using the following 
equations: 



10 

and 



where Y(x) is tiraa-diffcrentiated bioimpedance at a given point x, A, T, and B arc 
check point positions in each caxdiocycle (see FIG. 9); ST is the time between 

15 pointii S and T; and i and j are different cardiocyclcs. The computer system 
compares all of ihe calculated nearness vectors and eliminates those cardiocydes 
wherein the amplitudes of the. nearness vectors exceeds certain thresholds. Tlie 
comparison is made using a two-threshold analysis. If the nearness vector 
amplitude of two cardiocydes is less than a first predefined value, LI, the similarity 

20 is considered "good." If the similarity fails to be considered "good" (the similarity 
exceeds LI), but is still less than second predefined value, L2, it is considered as 
'acceptable. " If the similarity fails to be considered "acceptable" for the id second 
block of data, the computer system compares the tested cardiocycle with np to 50 
previous "not very bad" cycles. If there is still no cycle sijnilar to the one being 

25 considered, the cycle under consideration is regarded as noisy, and is rejected. If 
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the anioiiiu. nf "good" cardiocycleii is big enough, aU "acceptable" cardiocycles also 
rejected from ihc final calculation. This methodology increases the suhiJity of the 
calculations at the highest noise levels, as only "good" cardiocycles arc used in the 
final calculations. 

5 After the bioimpedance signal processing is complete, the computer system 

perfomis the final calculation of hemodynamic parajneters, together with means and 
variances for tlic entire data block. Heart Rate (HR) and Hean Stroke Volume (SV) 
are recalculated using the respective methods described above, except that the 
processed bioimpedance signal is used instead of the unprocessed bioimpedance 
10 signal- Cardiac Output (CO) is calculated as the product of HR and SV; that is: 

co^srsR 

An outline of the general methodology of the present invention is set fonh in 
FIG. 6 A. FIG. 6B depicts the flow chart of signal processing in the present 
invention. 

^5 While the present invention has been described in terms of a preferred 

embodiment, those of ordijoaxy skill in the ait will recognize and appreciate that it is 
not so limited. Many additions, deletions and modifications to Che disclosed 
embodiment may be made without dqjaiting from the scope of the invention as 
hereinafter claimed. 
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CLAIMS 

What is claimc<t is : 

1. An apparatus for detennining a subject '5 heqrt rate, heart stroke 
volume, and cardiac output from detected thoracic bioimpedance signals and 

5 electrocardiognun comprising on electrode array for detecting said thotacic 
bioimpedance signals and said electrocanliogTam. 

2. I'lie apparatus of claim 1 further comprising: 
an upper influencing electrode placed on the subject's head; 

10 a lower influencing electrode placed on the left lower extremity of the subject; 
an upper pair of detecting electrodes placed on the subject's neck; and 
a lower pair of detecLing electrodes placed on the trunk of the subject. 

3. The apparatus of cLaim 2, whcrcin the placement geometry Anther 
15 comprises: 

an upppjr influencing electrode placed on the subject's forehead; 
a lower influencing electrode placed in the general area of the subject's left knee; 
a pair of upper detecting electrodes placed on the subject's neck; and 
a pair of lower detecting electrodes placed laterally on opposite sides of the 
20 subject's chest. 

4. Tlie apparatus of claim 3, wherein said upper influencing electrode 
comprises a spot electrode for orientation on vertical and horizontal center lines of 
said subject's forehead. 



25 



30 



5. The apparatus of claim 3, wheniin said louder influencing electrode 
comprises a spot electrode, the placement of which satisfies the relationship L < 
5R, where L Is the vertical distance between said upper and said lower influencing 
electrodes and R is the XBdius of said subject's chest. 

6. The apparatus of claim 3, wherein said upper detecting electrodes 
compri-ses a pair of si)ot electrodes oriented symmetrically on oppo-sile sides of said 
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subjccCs iicclc along a horizontal line approximately 4 centimeters above the base of 
said subject's neck. 

7. The apparatus of claim 3, wherein said lower detecting electrodes 
5 further compri.se a pair of electrode assemblies, each assembly providing contact 
surface area between about 12 square centimeters and about 30 square centimeters 
oriented laterally on opposite sides of said subject's chest at appioxiinately xiphoid 
process level. 

10 8- The apparatus of claim 7, wherein each assembly fimher comprises 

four spot electrodes, wherein each spot electrode comprises approximately 4 square 
centimeters contact surface area with each said spot electrode centered at comers of 
a square with sides measuring 5 ccnthncters, and all 4 spot electrodes of said 
assembly electrically connected to each other. 

15 

9. The apparatus of claim 8, wherein the top spot electrodes of each 
assembly lie on the xiphoid process level of the subject. 

10. A method for processing a hioimpedance signal and 

20 electrocardiogram for deriving heart rate, heart stroke volume, and cardiac output 
comprising; 

registering gain-phase-frequency (GPF) characteristics of input analog devices for 

measuring bioLmpedancc; 
registering gain-phase-frequency (GPF) chaiacteristics of input analog devices for 
25 ujcasuring electrocardiogram; 

measuring bioimpedance as a function of time over a given time period with said 

bioimpedance input analog devices and generating a bioimpedance signal; 
measuring electrocardiogram as a function of time over aoid given time period with 

said elcctnocardiogram input analog devices and generating an 
30 elecLTOcardiogram (ECG) signal; 

correcting the bioimpedance signal for distortions based on the GPF characteristics 

previously rcgistenid; 
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conrcctnifi the clectrocardiogmm signal for distortions based on the GPF 
characteristics previously registered; 

determining valid QRS complexes associated with each caidiocycie of said 
electrocardiogram signal in the given time period; 
5 locating check-points on said valid QRS complexes; 

processing one of said ECG signal and said corrected bioimpedance signal Co 
itstimarr heart rate; 

time-differentiniing the concctcd bioimpedance signal; 

determining chcck-poinis for each said cardiocycle of the time-differentiated 
10 bioimpednnce signal in the given time period; 

determining effective left ventricular ejection time (ELVET) using said time 

differentiated bioimpedance check-points in relation to said corre^onding 
QRS check-points; 

determining a novel correction factor Z.^ using said time differentiated 
15 bioimpedajice check-points in relation to said corresponding QRS check- 

points; 

calculating stroke volume as a function of said ELVET, maximum time- 

differentiated bioimpedance (dZ/dt)^, specific blood reaistivity (P), distance 
(L) between two bioimpedance voltage sensing electrodes of the 
20 bioimpedance analog input device^ baseline bioimpedance (Zc) , said 

correction factor Z,^, and a novel scale factor (K); and 
calculating carch'ac output by multiplying said stroke volume by said heart rate, 

11. The method of claim 10, therein said registering gain-phase- 
25 frequency (GPF) characteristics of said bioimpedance input analog devices 

comprises determining phase-frequency and gain- frequency characteristics of a 
transducer employed in detection of said bioimpedance prior to use thereof in said 
detection. 

30 12. l^hc method of claim 10. wherein said correcting of measured 

bioimpedance signal comprises digitally fjlterins and phase correcting said measured 
bioimpedance to remove distortion in the output of said transducer. 
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13. 'I'hc method of claim 10, wherein said registering gain-phase- 
frequency (GPF) chamcteriscics of said elect rocardio gram input analog devices 
comprises dctcnnining piia^c -frequency and gain-frequency chaTacteristic^ of a 
transducer employed in detection of said electrocardiograxn piiur to use ihersof in 
5 said detection. 



10 



14. Tlie method of claim 10, wherein said correcting of measured 
elcctrocardiognim signal comprises digitally filtering and phase correcting sair< 
measured electrocardiogram to remove distortion in the output of jaid transducer. 

15. The method of claim 10, wheneia said estimating heart rate comprises 
using a power .spectrum of the bioimpedance signal, and an auto-convolution 
function of the said power spectrum. 



IS 16. The method of claim 10, wherein said estimating heart rate comprises 

processing the electrocardiogram signal. 

17. The method of claim 12, wherein said correcting of measured 
bioimpedance signal further consists of suppressing breath waves to remove 
20 undesired power spectra components and generate a bioimpedance signal of restore 
shape. 



18. The method of claim 10, wherein said determining valid QRS 
complexes comprises determining a distribution of all peaks measured in an 

25 electrocardiogram over a period of time (from the front (E,) and back (EJ 

amplitudes of the peaks, calculating the (E,, Ej) amplimdc envelope and rejecting all 
pcaics outside ilie envelope. 

19. The method of claim 11, wherein said determining transducer phase- 
30 fn:queiicy and gain-frequency characteristics comprises: 

generating a high precision sinusoidal impedance signal with peak-tcpcak 
impedance of approximately 0.2 Ohms and baseline impeilance of 
approximately 100 Ohms to 2U0 Ohms; 



* 
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connecling said sinusoidal iinpedajjce signal to said Iransdiicer; 
measuring the. (uicput from said transducer; and 

calculatiufi a gain-phase-frequency characterlsiic, H(f), of the transducer in a 
prederiiif.d frequency range. 

5 

20. Tiie method of claim 19, EiJrther comprising generating said high 
precision sinusoidal impedance signal using a voltage-to-impedance converter 
including a photoiesistor, a photoemitter, a power supply and an analog-to^igital-lo 
analog computer inieiface. 
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21. The method of claim 20, further comprising producing through said 
interface a set of tast signals with predefined amplitude and frequency range, 
detectuig an output signal from said tran.sduccr. and analyzing said output signal to 
dciciminc said phase-frequency and said gain-frequency characteiistics. 



22. nie method of claim 19, further comprising employing posterior 
signal processing to coneci linear gain-phaae-frcquency distortions by converting 
real operating characteristics of the transducer to predefined chaiacterisiics, wbenain 
phase shift is zeroed and gain is assumed to be constant in a predefined frequency 
20 range. 



23. A method of heart rate estimation, comprising; 
calculation of a power spectrum of a bioirapedance signal; 
multiplication of said power spectrum by a selected amplitude-frequency function 
25 differentiate the signal and suppress breath hannonics; 

autoconvoluting the resulting power spectrum according to the formula 



ASm=PSa(/)PSaCif>PSa(3/)... ; 



30 



and 

determining a maximum amplitude value of autoconvoliition in a predefmed 
frequency range as an estimation of heart rate. 
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24. A inediod for dcicrmining cardiocyclci, comprising: 
filteiinK a bioinipcdancc signal to emphasize fronts of cardiocycles; 
calculating a tinie-amplitude envelope of said caidiocycles by analyzing the first five 

hannonics of ihe power specLmm of said bioimpedance signal after said 
5 filtration; 

selecting said aitdiocycle fronts by comparison with said calculated time-amplitude 

envelope; and 
rejecting erroneously-detected fronts. 



10 25. A method of selecting valid cardiocycles from coixected 

bioimpedance signals to eliminate cardiocycles having interference artifacts, 
comprising; 

detecting time iind amplitude relations referencing check points within individuals of 
a plurality of cardiocycles; 
15 comparing said time and amplitude relations between individuals of a said plurality 

of cardiocycles; and 
further examining selected cardiocycles which exhibit the presence of artifacts 
according to a plurality of comparison criteria. 

20 26, The method of claim 25. further comprising: 

constnicting a multi-dimensional vector for each selected cardiocycle; 
comparing said inulti-dimcnsional vector with such vectors for other cardiocycles 
and; 

rejecting the cardiocycles with vectors having no neighboring vectors when 
25 compan:d to last 50 valid cardiocycles and other candidate cardiocycles. 

27. A method of deriving effective left ventricular ejection time from 
measured biouiipedance signal and measured electrocardiogram signal, comprising: 
filtcrinjj said measured bioimpedance signal and suppressing breath waves therein; 
30 filtering said measured electrocardiogram signal; 
detecting a valid cardiocycle; 

calculaiing the lime-derivative of said bioimpedance signal Y(x); 
deiemining the maximum value of the time-derivative (dZ/dt),„„; 
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dettrmining cfltctivc ejection start time (S-poiri); 
determining effective ejection end time (T-poiiit); and 
calculating effecuve left ventricular ejection time (BLVET) as change in time 
between effective ejection sLaurt nme and end lime. 

S 

28. The method of claim 27, wherein determining effective ejection start 
time comprises: 

determining Ihe global maximum for a given valid cardiocycle of a time- 
differentiated bioimpedance signal and designating said maximum as point A; 
10 tracing back in lime from corresponding point A on electnxardiogiara to point S,; 
looking for abuormalitles in the bioimpedance signal between points A and S^^; 
if there are no said abnormalities, then the cardiocycle is rejected as noisy; 
if there are any said abnormalities the one closest to point approached from the 
right is selected as the ejection start time S; 
15 otherwise the abnormality nearest approached from the left is selected as the 

ejection start tinic S. 



29. lite method of claim 28, wherein determining effective ejection end 
time comprises: 

20 determining the first (T,) and second (T,) local minimums at the time-differentiated 

bioimpKilance signal after point A; 
analyzing the depth of the signal curve ac each of the first (T,) and .second (TJ local 
minimums; 

if the depth of the second (T,) minimum is greater than a predetermined fractional 
25 value of the depth of the first (TO minimum, selecting the second minimum 

(T,) as To; 

otherwise, selecting T, as Tq; and 
identifying tlic T-point as the nearest local minimum before point Tq on the graph of 
the curve generated by the second derivative of Y(x). 

30 

30. The method of claim 29, further comprising, after idemifying T, and 
Tj and before identifying the T-point and regardless of relative ampHuidcs of Tj and 
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detectiiig a back edge of a T-wave in the electron cardiogram signal; 
determining if T, or Tj is out of bounds of the back of ihc T-wavc; and 
if so. .selcctint I he- one of Tj or which is not out of bounds 2S Tq. 

5 31. in the mctliod of claim 28, abnormalities are selected from the group 

compri.sing: d7./dt zero crossing (Q-poinO) local miiuiiium in dZ/dt| and local 
maximum in the third lime-derivative of the bioimpedance signal, d^Z/dt^. 

32. A method of determining stroke volume for a patient, comprising: 
10 determining specific blood resistivity P; 

measuring a distance L between two bioimpedance electrodes applied to the patient; 
determining tlic base thoracic impedance Zq; 
determining EL VET; 

determining AZ, impedance changes due to blood influx; 
15 and calculating stroke volume SV according to the equation 



where K is a novel scale factor related to body composidon of the patient. 

20 

33. The method of claim 32, further comprising calculating K as 

K = - • {SCHESTKH\ - W^)\ 

where 

SCHJEST = (PCHEST^^PNECKPCffEST^PNlECK^yilT^. 

25 34. Tl\e muihod of claim 33, wherein Kq, K,. Kj, K3 arc gender and age 

dependent and lie in ranges of 
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€ [1-4]; K, e [3-16]; X, e [0-11; e [0.1-2]. 

35. A method of daim 32, wherein deierrnining aZ comprises: 
locating the sun of QRS complex of an ECG signal and labeling it point Q; 
decennining the impedrmce ai point S, Z^; 
detennining the impedance at point Q, 21,; 
5 calculating the impedance difference Z,., between points S and Q; 
estimating AZ according to the formula 



AZ = idZJ df^^ ' ELVET ^ Z,,^ 

36. A method of breath wave suppression for a bioimpedancc signal, 
comprising: 

10 calculating the Fourier ccansforai of the signal; 

locating the first and second frequency harmonics of cardiocycles in the calculated 

spectrum of the signal; 
estimating the width of each of the harmonics; 

suppressing frequency harmonics below the lower bound of the second hanxionic 
15 except for harmonics within the bounds of the first frequency harmonic; and 

calculating the inverted Fourier transform of the signal. 
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